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I. INTRODUCTION 



In several recent papers [1-4] we have developed a highly predictive model of quark and 
lepton masses based on the grand unified group 50(10). This model grew out of our attempt 
IjU to construct a realistic grand unified theory (GUT) in which 50(10) was broken down 
to the standard model in the simplest possible, or "minimal" way 0. In this model there 
emerged a new mechanism based on certain well-known features of SU(5) for explaining 
the large mixing between the mu and tau neutrinos that is seen at SuperKamiokande ||. 
In HH we gave the structure of the quark and lepton mass matrices for the second and 
third families, treating the first family as massless. In it was shown how to extend the 
model to include the first family, which leads to several interesting predictions. In [|j], it was 
observed that the mixing of the electron neutrino very naturally falls either within the range 
0.004 < sin 2 29 efl < 0.008, corresponding to the small angle MSW solution of the solar 
neutrino problem, or very near to the value sin 2 # e/J = 1, corresponding to what is called 
"bimaximal mixing" . 

In this paper we present the model in fuller detail, especially in regard to neutrino 
phenomenology, and to the structure of the Higgs sector, Yukawa interactions, and flavor 
symmetries that underlie the quark and lepton mass matrices. 

The paper is organized as follows. In Section II we discuss in general terms, that is 
apart from a particular model, our mechanism for explaining the large mixing of and 
v T . In Section III, we explain what we mean by a "minimal" 50(10) breaking scheme, and 
show how such minimal breaking and the requirements of simplicity lead one naturally to a 
certain form for the mass terms of the heavier two families of quarks and leptons. We then 
observe that this form realizes the general mechanism for large — v T mixing described in 
the previous Section. It is important to emphasize that this mechanism emerged not from 
an attempt to explain neutrino phenomenology, but from other considerations entirely, in 
particular the attempt to simplify the Higgs structure of 5O(10). It is most interesting 
that the same mechanism has also independently been found by other groups attempting to 
make sense of neutrino phenomenology. In Section IV, it will be explained how this model 
is best extended to the first family of quarks and leptons, and how this gives rise to several 
distinctive predictions. Accurate analytic expressions for the predictions at the GUT scale 
will be presented. In Section V, the neutrino sector will be examined in detail. It will be 
seen how either the small-angle MSW solution of the solar neutrino problem or bimaximal 
mixing can result with equal simplicity. Finally, in Section VI, a concrete model, including 
all the details of flavor symmetries and of the Higgs and Yukawa superpotentials, will be 
presented, showing that the basic scheme is technically natural. 

II. MECHANISM FOR LARGE - v T MIXING 

Before explaining our mechanism, it will be helpful to explain why the observed large 
mixing of z/^, presumably with z/ T , has been a theoretical puzzle. The basic reason is simple: 
the mixing that is seen between the quarks of the second and third families is described 
by a small mixing angle, namely V c b = 0.04, and therefore it was expected that the mixing 
between the second and third family of leptons would also be small. 
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The grounds for this expectation were twofold. First, there is the empirical fact that 
the masses of the quarks and leptons exhibit roughly similar "hierarchical" patterns, and 
therefore it was natural to assume that their mixing angles would be similar also. Second, 
the most promising theoretical approaches to understanding the pattern of quark and lepton 
masses, namely grand unification and flavor symmetry, tend to treat quarks and leptons in 
similar ways. For instance, small quark mixing angles might suggest an underlying funda- 
mental "family symmetry" that is weakly broken, in which case the lepton mixings might 
be expected also to be small. And in grand unification based on 5*0(10) there is a close 
connection between the quark and lepton mass matrices. 

There are actually two puzzles associated with the mixing of the second and third fam- 
ilies: First, why is the lepton mixing [U^l ~ 0.7 so large? And, second, why is the quark 
mixing V c b = 0.04 so small? What we mean by saying that these are distinct puzzles is that 
they are both unexpected within the most commonly assumed framework for explaining 
quark and lepton masses, the Weinberg- Wilczek-Zee-Fritzsch (WWZF) idea ||. 

The WWZF idea was that the Cabibbo angle could be understood if the mass matrices 
of the first and second families of quarks had the following form: 

L mass = fecH) (° *J (^J+(S,35)(J *,J (fy. (1) 

This gives \md/m 8 \ = \b'/a'\ 2 , \m u /m c \ = \b/a\ 2 , and V us = b' /a' — b/a, and thus the famous 
relation 

V us = \Jm d /m s - e %a \]m u /m c . (2) 

Since \V US \ = 0.22, ^Jmd/m s = 0.22, and \Jm u /m c = 0.07, this relation is satisfied for 
a ~ ±7r/2. If we apply the same idea to the leptons of the first two families we get 

U e2 = ^rrie/m^ - e %p ^m Ul /m V2 . (3) 

The second term on the right is not known, but if it is assumed to be small one has the 
rough prediction that U& ~ ^rrie/m^ = 0.07. This could be consistent with the small angle 
MSW solution of the solar neutrino problem, which requires that U e 2 ~ 0.04. Thus the 
WWZF idea appears to work well where it was originally applied, namely to the first and 
second families. 

Fritzsch || later extended this idea to explain the mixing of the third family. If a WWZF 
form is assumed to hold for the second and third family, i.e., if one takes (u, c) — > (c,t) 
and (d, s) — > (s, b) in Eq. (1), one obtains 

V cb = \jm s /m h - e ll \Jm c /m t (4) 

and 

= \frr^Jrn T - e l5 ^Jm U2 /m U3 . (5) 

Since ^Jm s /m h = 0.14, and J~mjm t = 0.04, one sees that the observed value of V c b = 0.04 
is too small by a factor of three or so. Assuming that the neutrino mass ratio in Eq. 
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(5) is small, and given that ^Jm^/rrir = 0.24, one sees that the nearly maximal value of 

~ l/y/2 = 0.7 that is observed is too large by a factor of three or so. 

These Eqs. (2-5) are based on the assumption of a hierarchical and symmetric form 
for the mass matrices. A key feature in our mechanism for understanding the large mixing 
of the tau neutrino is that it involves highly asymmetric mass matrices. As we shall see, 
the assumption of asymmetric mass matrices naturally explains why is larger than the 
Fritzsch value and V c b is smaller than the Fritzsch value by approximately the same factor. 

Consider, a toy model with SU (5) symmetry, which has a set of Yukawa terms of the 
following form: A33(5 3 10 3 )5h + A 2 3(5 2 10 3 )5h + A3 2 (5 3 10 2 )5h, with A 32 < A 23 ~ A 33 and the 
subscript H denoting a Higgs representation. These terms yield the following mass matrices 
for the second and third families of down quarks and charged leptons: 

^)m d + M (^)(^)m d , (6) 

with e < (7 ~ 1. Here we have labelled the fermions with a family index, instead of 
the names s, b, fi, and r, since the mass matrices in this case are far from diagonal. A 
crucial point to notice is that the matrix for the leptons, which we will denote by L, is the 
transpose of the matrix for the down quarks, which we will denote by D. This is a feature 
of minimal SU(5). It arises from the fact that the 5 representation of fermions contains 
the left-handed leptons, II, and the charge conjugate of the right-handed down-quarks, d R , 
while the 10 representation of fermions contains the charge conjugate of the right-handed 
leptons, In, and the left-handed down-quarks, d^. Thus, SU(5) relates D to L, but only up 
to a left-right transposition: D = L T . 

The transposition feature of SU(5) unification appearing in Eq. (6) results in the large 
element, a, of L producing an 0(1) mixing of Z 2 l with for the leptons, while in D for 
the quarks it produces a large mixing of the right-handed fields <i 2 i? and d^R- The mismatch 
between the large 1 2 l — mixing and the v 2 l — ^zl mixing, which is small (as will soon be 
seen), leads to a large mixing element. But the right-handed mixings of the quarks are 
not observable through standard model physics. What matters is the left-handed mixing of 
<i 2 L with d^L, which contributes to V c b, and is controlled by the small parameter e. 

The common statement that grand unification relates quark and lepton mixing angles, 
and thus V c b to U^, is very misleading. What is really true in general is that grand unifi- 
cation relates the mixing of quarks of one handedness to the mixing of leptons of the other 
handedness. Thus V c b and need not be directly related to each other. Of course, if 
the mass matrices are symmetric, as has almost always been assumed, the left-handed and 
right-handed mixings are the same, and hence V c b is directly related to U^. The most nat- 
ural interpretation, then, of the experimental discovery that \U^\ ^> |T4&| is that the mass 
matrices are highly asymmetric. This is the essential point first made in [10]. 

Not only does a highly asymmetric, or, as we will call it, "lopsided," form of the mass 
matrices explain the difference between the size of and V c b, but it also explains the fact, 
noted above, that is larger than the Fritzsch value and V c b is smaller than the Fritzsch 
value by about the same factor. The point is that the product of the two off-diagonal 
elements, e and a, is controlled by the fermion mass ratio. As is evident from Eq. (6), 
m s / m b — tt-2 ~ ecr. That means that the Fritzsch prediction for the mixing of di and Sl, 
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which is ym a /mb, goes approximately as y/ea. That shows that the Fritzsch prediction for 
the mixing angles is roughly the geometric mean between the true value of ~ cr and 
the true value of V c b ~ e. In other words, in our hypothesis of lopsided mass matrices, the 
surprising largeness of and the surprising smallness of V c b are two sides of the same coin. 

Another important feature of this mechanism should be emphasized. Almost all pub- 
lished explanations of the largeness of the — v T mixing trace it to some special feature or 
form of the neutrino mass matrix. Perhaps this is due to the purely linguistic fact that we 
talk about "neutrino mixing angles". But they could just as well be called the "charged- 
lepton mixing angles". They are really the angles expressing the mismatch between the 
neutrino mass eigenstates and the charged lepton mass eigenstates, just as the CKM angles 
are the mismatch between the up and down quark eigenstates. In our mechanism, the large 
value of is traceable to a peculiarity of the charged lepton mass matrix L, namely, having 
a large off-diagonal entry a. As we shall see in the next Section, having such a large entry 
helps to explain several other features of the quark and lepton mass spectrum. 

To sum up, the mechanism for explaining large v^ — Vr mixing proposed in |l]-|3| has three 
salient features: (1) the largeness of this mixing is due to the charged lepton mass matrix, 
which is (2) highly asymmetric, and which is (3) related to the transpose of the down quark 
mass matrix by 5f/(5). 

In the next Section we will see how a model with precisely these features arises very 
naturally in 50 (10) from very different considerations. 

III. FERMION MASS MATRICES IN MINIMAL SCHEMES OF 50(10) 

BREAKING 

The model that we shall examine in this paper emerged originally from our attempt to 
construct a realistic model based on 5*0(10) JTl[ in which 50(10) is broken to the standard 
model group, Gsm = SU(3) C x SU{2) L x U(l)y in the simplest possible way. We shall 
therefore start by explaining what we mean by minimal 5O(10) breaking. 

Since 5O(10) is a rank 5 group, it requires for its breakdown to Gsm at least two Higgs 
fields. One Higgs field is needed to break the rank of the group to 4, but this generally 
leaves an unbroken 577(5). The second Higgs field is needed to break 577(5) down to Gsm- 
The two breakings can occur in either sequence depending upon which Higgs field has the 
larger VEV and effects the first breaking. 

Whatever Higgs field gives superlarge mass to the right-handed neutrinos, as required 
for the standard seesaw explanation of the lightness of the left-handed neutrinos, will also 
break 5O(10) down to SU(5), and thus the rank to 4. There are two simple choices for this 
Higgs field: either an antisymmetric five-index tensor 126 or a spinor 16. In either case, one 
also expects a Higgs field in the conjugate representation, 126 or 16, to go along with it. A 
nice feature of the 126 is that this tensorial representation leaves unbroken a Z2 subgroup 
of the center of 5O(10) that acts as an automatic matter parity, whereas if a spinor Higgs is 
introduced, then matter parity is not automatic. On the other hand, to introduce 126 + 126 
is to introduce quite large representations that tend to make the unified gauge coupling go 
non-perturbative below the Planck scale, and that may be hard to obtain from superstring 
theory. In any event, it would seem that the use of a spinor-antispinor pair, 16 + 16, is more 
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economical. Thus we assume that the rank of SO (10) is broken at the unification scale and 
the right-handed neutrinos get mass from one such spinor-antispinor pair of Higgs fields. 

To break the group the rest of the way to Gsm requires the existence of Higgs fields in 
the adjoint representation 45 and/or in the symmetric two-index tensor representation 54. 
Most published realistic SO (10) models have several of both kinds of multiplets. However, 
it has been shown that it is possible to break SO (10) to Gsm with only a single adjoint 
Higgs and no larger representations |BJ. 

This, then, is what we call the "minimal breaking scheme for 50(10)": The breaking of 
50(10) to Gsm is accomplished by the expectation values of a set of Higgs fields consisting 
of 45 h + 16j? + 16#, with the model containing no multiplets larger than the single 45. 
There, of course, have to be other Higgs fields to break the SU(2)l x £7(1) y group of the 
electroweak interactions. 

This minimality assumption is restrictive enough that it is possible to say in which 
direction the expectation values of these fields point. This can be done by considering the 
problem of doublet-triplet splitting, whereby the colored partners of the weak-doublet Higgs 
fields of the standard model become superheavy while the weak-scale masses of the doublets 
themselves are preserved. In 50(10) the only known way of doing this in a technically 



natural manner is the Dimopoulos-Wilczek or "missing VEV" mechanism [|T^J. The idea 
is that if an adjoint Higgs field that has an expectation value proportional to the 50(10) 
generator B — L couples to Higgs fields in the vector representation, it will make their color- 
triplet components heavy (since they have B — L = ±2/3) while leaving their weak-doublet 
components massless (since they have B — L = 0). The needed coupling is simply of the form 
10ih45h102h- Of course, the expectation value of the adjoint, by virtue of the definition 
of the adjoint representation, is necessarily a linear combination of generators of the group. 
B — L is one of the 50(10) generators that is picked out by simple forms of the Higgs 
superpotential for the adjoint multiplet. It should be noted that there is another version 
of the missing VEV mechanism that works in 50(10), in which the VEV of the adjoint is 
proportional to the generator I 3R of the SU(2) R subgroup of 50(10) |fT3"| . However, that 



version is significantly more complicated. Therefore, simplicity dictates the choice that 

(45 H ) cx B - L. (7) 

Since the assumption of a minimal 50(10) breaking scheme included the supposition that 
only one adjoint exists in the model, no adjoint exists except the one that points in the 
B — L direction. As we shall see, this puts an important limitation on the possibilities 
for constructing realistic mass matrices for the quarks and leptons. The assumption of a 
minimal 50(10) breaking scheme thus acts as an important guide in searching for good 
models. 

The simplest possible terms that would give mass to quarks and leptons in 50(10) would 
be \ijl6il6jl0n, where the subscripts i and j are family indices. This would lead to four 
proportional Dirac mass matrices for the up-quarks (U), down quarks (D), charged leptons 
(L), and neutrinos (N). In fact one would have D = L oc U = N. Moreover, all these 
matrices would be symmetric, which is why one can write D = L instead of D — L T as in 
minimal SU(5). Some of the predictions that follow from these relations are good, notably 
the famous prediction = ml where the superscript zero stands for quantities evaluated 
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at the unification scale Mq. However, D = L also predicts that mP s = and m d = m° e . 
Empirically, one finds instead that m° s ~ and m d ~ 3m". These factors of three are 
called the Georgi-Jarlskog factors fI3 |. The simplest possible SO(10) Yukawa terms also 



predict that all the CKM angles vanish, since U oc D. While not exactly true, this is at 
least a good zeroth order relation, since the CKM angles are all small compared to unity. 
By contrast, in SU(5) the matrices D and U are not related by the unified symmetry and 
so the CKM angles are unconstrained. The smallness of the CKM angles can be regarded, 
therefore, as evidence for 5*0(10). On the other hand, the proportionality of D and U in 
50(10) also predicts that m° c /m® = m°/m°, which fails badly by over an order of magnitude. 

What one can conclude is that a way of going beyond the simplest possible 50(10) 
Yukawa scheme must be found which preserves some of its predictions while breaking others. 
One way to do this involves using larger representations to break the electroweak interactions. 
For instance, in the original Georgi-Jarlskog model, a 45 multiplet of SU (5) (not to be 
confused with the adjoint of 5O(10)) participates in breaking SU(2)l x U(1)y- In the 
context of 5O(10), this 45 is contained in a 126, which is inconsistent with our minimality 
assumptions. More economical is to assume that the Higgs fields that break 5O(10) at the 
unification scale, i.e., the 45// + 16// + 16//, couple to quarks and leptons and thus introduce 
the effects of that 50(10) breaking into the quark and lepton mass relations. This is the 
assumption we make. 

To describe the third family it is simplest to assume the minimal Yukawa term I63I63IO// 
as pictured in Fig. 1(a). By itself, this would make all the mass matrices have the form 



(8) 



That would give the following predictions, all of which are at least good zeroth approxima- 
tions to reality: = mj, V c b = 0, and m\jm^ = m2/m 3 = 0, where m 8 is a mass of a 
fermion of the i th family. Note that these are just the "good" 5O(10) predictions mentioned 
above. 

The second family presents more of a challenge. The main issue is how to get the Georgi- 
Jarlskog factor of 3 between and mP s . Breaking of 577(5) must be involved, since the bad 
relation m° s = m° arises already at the 577(5) level. The only field that breaks SU(5) in the 
framework of minimal 5O(10) breaking is the adjoint, 45//. Since (45//) oc 77 — 7, and the 
77 — L of leptons is —3 times that of quarks, this field has the possibility of giving the needed 
Georgi-Jarlskog factor. Thus one must seek an effective Yukawa term that involves the 45//. 
The simplest such term ||15|| , in the sense of the term of lowest dimension, is of the form 
(16jl6j)10//45///MG'. Moreover, this term can arise in a simple way by the integration out 
of a 16 + 16 family plus antifamily at the unification scale, as shown in Fig. 1(b). 

There are actually two ways to contract the 5O(10) indices of such a term: the product 
(16jl6j) can be contracted symmetrically or antisymmetrically. It is easy to show that if 
(45) oc 77 — L, only the antisymmetric contraction contributes to the quark and lepton mass 
matrices. (The reason is simple. If the VEV of the adjoint is proportional to a generator 
Q, then the symmetric/antisymmetric contractions give contributions to fermion masses 
that go as Q(f) ± Q{f)- Since B — L of an antifermion is minus that of the fermion, the 
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contribution cancels for the symmetric contraction.) Thus, one need only consider the flavor- 
antisymmetric term, which means only ij = 23 and not ij = 22, or 33. Consequently, the 
only operator of interest is (16 2 16 3 )10^45^ which, together with the operator I63I63IO//, 
gives 



U 



N 



(0 



















6/3 










( : 





— e 


\o 


e 


1 







D 



M, 



u, 



(0 











Vo - 


-6/3 


fO 








— e 


VO e 


1 







(9) 



ID- 



The desired factor of 3 has been achieved between leptons and quarks, due to the gen- 
erator B — L to which the adjoint VEV is proportional. One also can see that the e entries 
are flavor antisymmetric for reasons already explained. As they stand, these forms of the 
matrices are inadequate to explain even the features of the second and third families of 
fermions. There are three inadequacies. (1) The factor of 3 comes in squared between the 
mass of the leptons and quarks of the second family. The reason is that, for e small due 
to the mass hierarchy between families, the second eigenvalue of L is given by the seesaw 
formula m° = e 2 Mo, while the second eigenvalue of D is given by m' = (e/3) 2 M£). (2) The 
matrices D and U are still exactly proportional. This is a consequence of the fact that the 
generator B — L does not distinguish up and down quarks. Therefore, the CKM angle V c b 
still exactly vanishes. (3) Because D and U are exactly proportional, one still has the bad 
prediction mj)/m t ° = m°/m°. 

It is clear that the breaking of SO (10) due to the adjoint cannot cure all of these problems, 
since B — L does not distinguish D from U. Thus the breaking of 5*0(10) done by 16# + 16# 
must come into play. As we shall now show, a single simple operator exists, which involves 
one of these spinor Higgs and cures at one stroke all three of the problems we have identified. 

The lowest dimension effective Yukawa operators that involve the spinor Higgs fields are 
quartic in spinors. Consider, therefore, operators of the form 16^16^ 16'^ h^h/Mq. The 
16// is the spinor Higgs field that breaks SO (10) at Mq down to SU(5). The 16'h is a 
spinor Higgs that has a weak-scale VEV that breaks SU(2)l x U(1)y- In principle, these 
two spinors could be the same field. However, if they were, it would mean that they had 
to be contracted symmetrically by Bose statistics, which in turn would mean that 16j and 
16 j would also have to be contracted symmetrically. A careful examination shows that 
the resulting flavor-symmetric contributions to the mass matrices do not lead to realistic 
forms, though it is possible to achieve realistic mass matrices by adding yet another Yukawa 



operator, as in the interesting model of Babu, Pati, and Wilczek [fT6 ]. Therefore 16' H must 
be a distinct field. As will be seen later, introducing this 16' H involves no loss of economy, 
since it allows a very elegant explanation of the largeness of the ratio m t /mi, without making 
tan/3 large. 

There are still several operators of this type to be considered: the family indices can 
take the values ij = 33, 22, 23, or 32, and there are three ways to contract the four spinors 
to make an SO (10) singlet. Here again, one must examine the various cases to see which 
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gives the most realistic mass matrices. As it turns out, there is one operator that is much 
superior to the others, in the sense that it much more cleanly and simply fits the data. It is 
of the form [I62I6//] [163l6'#], where [...] means that the spinors inside are contracted into 
a 10. This can arise very simply by integrating out a 10 of fermions, as shown in Fig. 1(c). 

Let us write the resulting mass operator in SU (5) language. Denote by p(q) a p multiplet 
of SU(5) that is contained in a q multiplet of S'O(IO). The VEV of 16h lies, of course, in the 
1(16) direction, while the VEV of 16' H that breaks the weak interactions lies in the 5(16) 
direction. Therefore, the resulting mass term is of the form [5(162)1(16#)][10(163)5(16 / // )], 
which in SU(5) terms gives effectively the operator (5 2 103)5h- Note that this has the same 
form as the SU(5) operator discussed in the last Section, which gave the a entries in Eq. 
(6). The result, then, of including this operator || is to make the mass matrices take the 
form: 
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The new term has given the entries we call a. Note that these lopsided entries appear only 
in D and L. The reason is simply that the 16' H contains a 5 of SU(5) but no 5. 

It is easy to see that the new term with o ^> e cures at once all three of the problems we 
identified with the forms given in Eq. (9): (1) Instead of mJJ = e 2 Mp and m° s = (e/3) 2 Mc, 
one has approximately that m° oc (e)(cr + e) = ecr and m° s oc (e/3) (a + e/3) = ea/3. More 
exact expressions will be given later. Thus the desired Georgi-Jarlskog factor of 1/3 is 
obtained, instead of 1/9. The a entry has dominated over one of the factors of e/3 and thus 
prevented the factor of 1/3 from coming in squared. 

(2) The a entry comes into D but not U, and thus breaks the proportionality of the two 

2 

matrices. As a result, V c b no longer vanishes, but is given approximately by (e/3)(^^qr[). 

Note that this is of the same order in e as m s /m b = (e/3)(^-j-), rather than ^Jm s /m b as is 
the case with Fritzsch forms, and accords much better with the actual experimental values. 

(3) The fact that a breaks the proportionality of U and D also means that the bad 
relation m°/m° = m°/m° is broken. Specifically, mPjm® is of order e, while mPjm® is still of 
order e 2 and therefore much smaller. This also accords well with the experimental numbers. 
In fact, as we shall see, if one uses V cb and m M /m r to fix the parameters a and e, one finds 
that m c (l GeV) is predicted to be in agreement with the experimentally determined value 
of 1.27 ± 0.1 GeV. It should also be noted that the prediction m° = mj is only very slightly 
affected by the addition of the a term, both m' and mf being given to leading order in e by 

The economy of the above mass matrix forms is seen in the fact that five quantities (V^, 
m fi/ m T, ™ s l m bi ™ c / m t-i and m T /m b ) are successfully fit with only the two parameters a and 
e. No other published form succeeds in accurately reproducing the masses and mixing of the 
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the heavier two families with so few parameters. The predictions and fits will be discussed 
in detail in Section IV. 

We see that the matrices in Eq. (10) were arrived at by a process of reasoning that had 
nothing to do with the question of neutrino mass but rather with an attempt to get realistic 
masses and mixings for the quarks and charged leptons using as simple a Higgs sector as 
possible in S'O(IO). But what has emerged is a structure with precisely the three critical 
features identified in the last Section as giving a simple explanation of the large mixing of 
and v T . In fact, a fit of V c b and m M /m T gives a = 1.8 and e = 0.14. Consequently, as can be 
seen directly from Eq. (10), the angle 6^ T = sin -1 = tan -1 a — 0(e) = 60° — 0(8°). This 
is quite consistent with what is observed. We will look more carefully at these predictions 
later. 

To summarize, with two parameters, e and a, four mass ratios and two mixing angles 
are satisfactorily accounted for, if we include U^. No greater economy could be hoped 
for in explaining the spectrum of the heavy two families. Moreover, as we shall see in the 
next Section, the forms in Eq. (10) can be extended to include the first family with equal 
economy: the introduction of two new parameters (one of which is complex) will nicely 
account for seven quantities pertaining to the first family. 

Before explaining how the model is extended to the first family we will expand on a 
couple of points made earlier. First, we said that the introduction of the 16' H allows a 
simple explanation 0,0 of why m t ^> m& that does not require a tan/5 ^> 1. The point is 



that the Higgs doublet of the MSSM that is often called Hjj is purely contained in the 10# 
that couples to I63I63 and gives rise to the "1" entry in the mass matrices of Eq. (10). 
However, the Higgs doublet of the MSSM that is called Hd does not come purely from 10#. 
Rather it is a mixture of doublets in 10# and 16' H , since they both contain 5's of 577(5). 
Thus we may write 

77^ = 77(10^), 

(11) 

H D = H(10 H ) cos 7 + H(16' H ) sin 7, 

where 7 is some mixing angle that depends on the parameters of the Higgs sector. Since 
the 33 elements of the mass matrices all arise purely from the coupling of the 10 h, the 
parameters we called M v and M D in Eq. (10) are given by 

M u = X 33 (Hu), 

(12) 

M D = X 3 3(H D ) cos 7. 

leading to the ratio Mjj/M d = {Hu) /{{Hd) cos 7) = tan/?/ cos 7, where tan/3 is defined to 
be the ratio of the Higgs VEV {Hu) giving mass to the top quark to the Higgs VEV {Hd) 
giving mass to the bottom quark. Hence 

m°/m° ^ (a 2 + l)- 1 / 2 (tan/3/ cos 7) (13) 

The point is simply that the large ratio of the top to bottom masses could be the result of 
cos 7 being small rather than tan/? being large. In fact, since we do not know anything a 
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priori about the angle 7, we cannot say whether tan (3 is large or small. It should be noted 
that if one assumes that Hp lies mostly in the 16'^ (so that cos 7 <C 1), it would explain 
why the parameter a is large since it comes from a coupling to 16' H , and also explain why 
tan (3 might be small. 

A second point we wish to underline here has to do with the reasonableness of asymmetric 
mass matrices. In many models it is assumed that all the mass matrices are symmetric. 
However, this is not something that is called for by the group theory of grand unification. 
It is true that with the minimal Yukawa terms SU (5) gives a symmetric U. But even with 
minimal Yukawa terms SU (5) does not predict any symmetry of the D and L matrices. 
And in 5*0(10), as we have seen, once one introduces the effects of 5*0(10) breaking into 
the Yukawa sector, as one virtually must, one easily obtains effective Yukawa terms that are 
asymmetric. Fig. 1(c) shows that very simple diagrams can give terms that are lopsided, 
in the sense that they contribute only above or below the diagonal. From the point of view 
of the fundamental grand unified theory, then, lopsided terms are as natural as symmetric 
ones. The preference for symmetric terms has been the result not of examining what kinds 
of terms are obtained in a simple way in unification, but rather from the desire to reduce the 
number of parameters at the level of mass matrices with the aim of making models which 
are highly predictive. However, putting oneself in the straight jacket of symmetric matrices 
makes it hard to get a good fit to all the quark and lepton masses and mixings. It turns 
out, as we have seen, and will see further below, that allowing asymmetric matrices makes 
possible a model which gives both a very good fit to the data and is actually much more 
predictive than most models which assume symmetric matrices. 



In arriving at the form of the mass matrices for the heavy two families we were limited in 
the choices that were possible by the assumption we made about the simplicity of the 50(10)- 
breaking sector. In extending to the first family we are not quite so limited. Nevertheless, 
the number of simple possibilities is not very large. There are several discrete choices: Should 
the contributions to the first row and column of the mass matrices be flavor symmetric like 
the l's in Eq. 10, antisymmetric like the e's, or lopsided like the ex's? Should they contribute 
to all the matrices equally like the l's, to all the matrices but with non-trivial Clebsch factors 
like the e's or only to D and L like the ex's? It is fairly easy to run through the various 
cases and see what kinds of relations among masses and mixings result. As it turns out, one 
of the simplest possibilities gives a remarkably good fit to the data. This uniquely simple 
choice |3| is the following: 



IV. EXTENSION TO THE FIRST FAMILY 



u 




5 5' \ 
a + e/3 M D , 
e/3 1 J 



(14) 
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We have already mentioned that fits give a = 1.8 and e = 0.14. The new parameters 5 
and 8' both have magnitude of about 0.008. The parameter 77 is by far the smallest, being 
about 8 x 10 -6 . The only role that 77 plays in the sector of quarks and charged leptons is 
in giving the up quark a mass, for it makes negligible contributions to the down quark and 
electron masses as determined from D and L, respectively. In Fig. 2(a) we have illustrated 
a higher-order diagram that can contribute to the parameter 77. Since it is not excluded that 
the up quark is exactly massless, it is possible to set 77 to zero. In any event, one can see 
that 77 = m°/mf, which is by orders of magnitude smaller than any other interfamily ratio of 
masses in the standard model. It will, however, be of some significance for neutrino masses. 
If 77 vanishes, this model gives only the small-angle MSW solution to the solar neutrino 
problem. But even if 77 is as small as 8 x 10~ 6 , it allows either the small-angle MSW solution 
or bimaximal neutrino mixing to arise in a simple way. 

Turning to the parameters S and 5', we see that they appear symmetrically and only 
in D and L. Such terms are easily obtained in SO (10) from simple diagrams such as 
that shown in Figs. 2(b) and 2(c). The effective operators arising from these diagrams 
are of the form [1 61 16j] [16^16'^] with j = 2,3, where again the spinors in brackets are 
contracted symmetrically into a 10 of 5*0(10) which is integrated out. Note, however, 
that the symmetric contributions 5 and 5' from the two Higgs contraction [1(16^)5(16^)] 
contributes only to D and L by virtue of their SU(5) structure. Contrast these effective 
operators for j = 2, 3 with that occurring previously for the term a arising from the diagram 
shown in Fig. 1(c). 

The three new parameters we have introduced are, as we shall see, sufficient to account 
for everything about the first family. Before proceeding, however, we must be careful about 
complex phases. It is easy to show that if we allow all parameters of the model to be 
complex all but two phase angles can be rotated away from the mass matrices U, D, L and 
N, provided we now neglect the negligible 77 contributions to D and L. We will call these 
physical phases a and (f> which appear as follows, 



U = 



N 




D 




( 





s 



S'e* \ 



(15) 



M 



D; 



) 



5 

-e 

y 5 y(0+ Q ) a + ee ia 1 

where in these matrices and henceforth e, a, 5, 5' and 7/ denote the magnitudes of these 
parameters and the phases are written explicitly. The phase a only comes into the fits of 
masses at higher order in the small quantity e/a = 0.08. Numerically, its effect is only a 
few percent; moreover, the fits (especially to m c ) prefer a value near zero. Therefore, we 
can ignore a and will do so from now on. That leaves only the phase 0. Its only significant 
effect, but a very important one, is to give the CP-violating phase angle 5cp- 

Instead of using the parameters 5, 5' and e 1 ^, it will be somewhat more convenient to 
use the parameters ix,, and e td , which are defined in terms of them as follows: 

, fl S-aS'e** . . 
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and 



The significance of these parameters is that they are essentially the left-handed and right- 
handed Cabbibo angles. This can be seen by taking the forms for D and L given in Eq. (15) 
and diagonalizing the 2-3 block. When this is done the 1-2 blocks of these matrices take the 
form 

^ <*(£?)• LM -(Vs)- < i8 > 

In terms of the five dimensionless parameters e, a, ti-, t R , and e l9 with r/ set equal to zero, 
we now write down expressions for fourteen observable quantities: seven ratios of quark and 
lepton masses, three CKM angles and one phase, and three lepton mixing angles. 



ml/ml - 1 - f^e, 
m°/m° = ±e 2 • [1 - fe 2 



9^ L 9 

mljml - • [1 + e^^f + ±{t\ + t R )\, 

<K = ■ [1 + | e ^rr^ + §(*£ + tDh 

m ° u / m t = 0, 

ml/ml - \t L t R ■ [1 - + e 2 ^±^ _ + 

= t L t R ■ [1 - |6^_ _ + t |) + (t 4 + ^ + t 4 }]) (1Q) 

V <*> — | e ^TT ' t 1 + l e a(J+l)]' 

1/ ~ f T [1 _ 1/2 _ / 2 , f 4 , 5 f 2 f 2 , 3,4 _ t R . -jfli 

v us — <-L[-L 2 i L i R ^r L R ^r 2 L L L R^ 8 L L 3aV^+l t L e J' 

C/°3 = sinV = ^T + C(e), 

^ 2 = cos V ■ [1 + e (^f - ^fl^Sfii) - if* _ l tll 

U e3 = tan 6^ ■ [1 + e -|^ fee*^ - 



These expansions have been carried to sufficiently high order in small quantities to be ac- 
curate to within 0.2% and are useful in doing the fits to the data. However, the leading 
terms in these expansions have much simpler forms and thus allow one to see more readily 
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the relationships among various quantities in this model. We therefore write these simpler 
expressions for purposes of discussion. 

m° b /m° T = 1, 
m° c /m° t = \e 2 



9^ ' 



,0 /™0 ~ , cr 
fc T 2 



a 2 + l- 



m° d /m° s t L t R , 



■ l f 



C = ^^(v^TTff - 1), 

C/°3 = Shi V = 7fe + 0(c)~ 0.7, 

C/ e ° 2 =008^(1^), 

C/ e3 = sin^ T (|^). 



(20) 



It might at first seem surprising that without any information about the Majorana mass 
matrix Mr of the right-handed neutrinos we are able to write down predictions for the 
three neutrino mixing angles. However, if rj = 0, as we are assuming at present, then the 
Dirac mass matrix of the neutrinos (N) has vanishing first row and column, and therefore, 
obviously, the same will be true of the mass matrix of the light neutrinos, which is given 
by the well-known "see-saw" formula M v = —N T M^ 1 N. This means that the two mixing 
elements of the electron neutrino, U e 2 and U e z, get no contribution from diagonalizing M u , 
but come entirely from diagonalizing L. Since L is a known matrix in our model, these two 
mixing elements are predicted. In the case of the mixing of the mu and tau neutrinos, 
does receive a contribution from diagonalizing M v . However, as can be seen from the form 
of N this is an effect of 0(e). The contribution to coming from diagonalizing L, on 
the other hand, is of order unity, since it arises from the large parameter o. Thus is 
predicted, although not precisely. 

Since we have written fourteen quantities in terms of five parameters, there are altogether 
nine predictions of the model. Which quantities one takes as "predicted" depends on which 
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quantities are used to determine the values of the parameters. We will use the lepton mass 
ratios and the angles V c b and V us for this purpose as they are the best measured. As one 
can see from the third and eighth of Eqs. (20), one can get the value a from the ratio 
31^°/(m°/m°). One finds (of course, taking the renormalization effects into account as was 
done in [g]) that numerically a ~ a/3. Substituting this into the expression for m"/mj, one 
obtains that e ~ 0.14. This is the small parameter of the model that is responsible for the 
hierarchy between the second and third families, and is small enough that the expressions in 
Eqs. (20) are fairly accurate. One can use m e /m^ and V us to determine ti and t R . A careful 
fit, described later, gives t L = 0.236 and t R = 0.205. That t L ~ t R is easily understood from 
Eq. (18) and the well-known Weinberg- Wilczek-Zee-Fritzsch result || that the Cabbibo 
angle is well acounted for by symmetric mass matrices for the first two families; cf. Eq. (1). 
The near equality of and t R is also apparant from the seventh and ninth relations of Eqs. 
(20) and the fact that numerically V us = \jrn^fm s . The phase factor e td will be determined 
from the CP- violating phase 5c p- 

The nine predictions, then, are the following. To begin with, there are the three famous 
predictions, (1) m° b /m° T = 1, (2) m° = §m°, and (3) m° d ^ 3m°. The first is the "good" 
prediction of minimal SU(5) unification, and the latter two are the Georgi-Jarlskog relations. 
These predictions are manifest from the first, third, fourth, sixth, and seventh of Eqs. (20). 
It is hardly surprising that the model gives these relations, since we were guided by them 
in constructing the model. The fourth prediction is (4) m°/m° = 0. Even if the u quark 
is not exactly massless this relation is a very good approximation to reality. If one takes 
the favored value of m u « 4 MeV, then, with reasonable assumptions about thresholds in 
doing the running up to the GUT scale, one obtains r] ~ mj/mj « 8 x 10~ 6 . This is far 
smaller than any other interfamily ratio of masses. For instance, the comparable ratio for 
down quarks is mj/mj = 10~ 3 , and for leptons is m° e /m° T = 3 x 10~ 4 . Like the previous 
three relations, m u « is a reflection of basic group-theoretical aspects of the model. It 
comes from the fact, explained above, that the 5 and 5' entries only appear in D and L. 

The remaining five predictions are not simple group-theoretical relations like the forego- 
ing, but are non-trivial quantitative predictions. They are predictions for (5) m c , (6) V u b, 
(7) C/ M3 , (8) U e2 , and (9) U e3 . 

The prediction for m c is particularly interesting. We see immediately that, for reasons 
having to do with the group-theoretic structure of the model, the ratio m° c /m® is much less 
than the corresponding ratio m°/m° for the down quarks because it is of higher order in the 
small parameter e. This is a highly significant success, because the minimal Yukawa terms 
of S'O(IO) notoriously give these ratios to be equal. Moreover, the success is not merely 
a qualitative one. When e and a are fit (using V c t and m^/m T ) and the renormalization 
effects are later taken into account, it is found that m c comes out within about 5% of the 
experimentally preferred value, which is quite remarkable given the various experimental 
and theoretical uncertainties. This success is non-trivial, because the reasoning that led to 
the forms of the mass matrices did not depend upon the value of m c , and hence it could 
have been expected that m c would come out wrong by a large factor. 

Another non-trivial quantitative hurdle for the model is the prediction for V^. The 
eighth, ninth, and tenth relations of Eqs. (20) give V° b = V^ s V^{y/a 2 + l^e~ i0 - 1). If 

we use the facts that a ~ \/3 and ti ~ t R , this gives V u b — V us V c b(^e~' t9 — |). A careful fit 
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gives 

V ub = KsK fe (0.558e^ e - 0.315). (21) 

In other words the model predicts that V u b should lie on a certain circle in the complex 
plane. As can be seen from Fig. 3, the circle for V u dV* b slices neatly through the middle of 
the presently allowed region. Again, this is a very significant success, since the reasoning 
that led to the forms in Eq. (15) was not based on the value of V u b- 

The prediction for the mixing of and v T has already been discussed. It is one of the 
key successes of this model that this mixing turns out to be nearly maximal. The fact that 
a ~ tells us that the first term in the expression for C/° 3 in Eq. (20) corresponds to an 
angle near 7r/3. As we shall see in the next Section, the 0(e) corrections easily bring this 
down close to the maximal mixing value of 7r/4. 

The prediction of this model for the mixing of v e and with rj = is quite interesting. 

From the sixth relation of Eqs. (20) and the fact that ti — £r, one sees that ~ ^Jrrie/m^. 
Thus the model predicts that U e 2 — cos 6^ T ^m e /m^. The factor of cos 6^ T is crucial fTF 



since without it one would have sin 2 26*^^ = 4|£/ e i| 2 |£7 e 2| 2 ~ 4(m e /m /J ) = 2 x 10 -2 , which 
is about twice the value needed for the small-angle MSW solution to the solar neutrino 
problem. Since atmospheric neutrino data tells us that cos 6^ ~ l/v^, the model gives just 
the correct value for the small-angle MSW solution. 

In the future both V u b and sin 2 26' so i ar will be known better and will provide a sharp test 
of the model. The theoretical uncertainties in the predictions for V u \, and U e 2 are estimated 
to be only a few percent. 

In discussing the v e — mixing above, we have assumed that 77 = 0. If 77 does not vanish, 
but is around 8 x 10~ 6 , corresponding to m u « 4.5 MeV, then it turns out that both the 
small-angle MSW solution we just discussed and bimaximal mixing are possible. This will 
be discussed in detail in the next Section. 

Finally, there is the prediction of v e — v T mixing. One sees from Eq. (20) that there 
is a prediction that f/ e3 = t&n9^ T U e 2 = 0.05. It is interesting that even for the bimaximal 
mixing case that will be discussed in the next Section, the numerical value of U e 3 is virtually 
unaffected. Thus this prediction is a "robust" one of this model. 



V. NEUTRINO MIXING 
A. Mixing of — u T 

In analyzing the predictions of this model for z/ M — v T mixing, we may make the approx- 
imation that 7] = 0. This means that N has vanishing first row and column. Therefore, in 
computing M u = —^M^N the first row and column of M^ 1 are irrelevant. Thus we may 
write M^ 1 as 



M~ R l = - X Y , (22) 
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where X, Y, and Z are in general complex. There are consequently five real parameters (the 
over all phase does not matter) that come into the masses and mixing of and u T from 
Mr. As observed earlier, this does not prevent us from making a qualitative prediction for 
the mixing parameter U^, since the contribution to it from diagonalizing the mass matrix 
M v is only of order e and comes predominantly from diagonalizing the known matrix L. 
However, in order to see if more precise predictions can be obtained, we shall look at two 
simple special cases: 



(/) M J 



R 



I - 

Be if3 e 
\ Be lf3 e 1 



A R e 



»7 



(23) 



and 



(//) M t 



R 




(24) 



In these cases only three parameters in Mr, namely Ar, B and e l @, contribute to the neutrino 
observables of the second and third families, since e has appeared previously and is used as 
a natural scaling parameter. In the first case, 



Mt = - 



/0 
e 

\0 e 2 + B- 1 e~ il3 



M 2 

m u c -i(0+~t)_ 



BA 



(25) 



The neutrino mixing matrix U, now known as the MNS mixing matrix fl8|| , is given by 
U = U\U V , where Ul is the unitary matrix that diagonalizes L^L, and U u is the unitary 
matrix that diagonalizes M\,M V . For case I, U v is given by 



10 

cos 6% sin ^ 3 
v -sin^* cos ^ 3 



(26) 



where tan 2 #23 



2e/K, and K = 2+B x e l/3 . The ratio of eigenvalues of M v gives m U2 /m 



(e 2 / \K\ )(1 — e 2 / \K\ + ...). One can choose cos #23 to be real, and one can write 



sin ^3 = 



m V2 /m v 



1 + 



/m 



(27) 



where e 1 ^ is the phase of K. One readily sees from the form of the charged-lepton mass 
matrix L in Eq. (15) that sin ^3 = (^l) 23 is given by tan2^ 3 = - J^X e = -J^ + 0(e). 
Since a ~ v3 it is evident that 6*33 — 60°. Using the best fit values of a and e one finds, 
more precisely, that #23 = 63°. 

Altogether, then, the mixing parameter of and v T is given by 
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U^ 3 = sin 6^ 



- Sin #23 cos $23 + cos $23 sm $23 



(28) 



= -0.898(1 - m v Jm V2 ) + O.Ul^m U3 /m U2 i 



If neutrino masses are hierarchical, and atmospheric neutrino oscillations are — v T oscil- 
lations, then m V3 ~ 0.06 eV. And if one further assumes the small-angle MSW solution to 
the solar neutrino problem, then m U2 ~ 0.003 eV. Thus, m U2 /m U3 ~ 0.05, within a factor 
of two or so. Taking it to have the value 0.05, and the phase £ to vanish, Eq. (28) gives 
U^ 3 = —0.756, and sin 2 2# MT = 0.984. With the same value of the neutrino mass ratio and 
£ taken to be 7r/4, sin 2 29^ T = 0.943. We see that there is excellent agreement with the 
experimental limits from SuperKamiokande if the complex phase is not too large. But if 
£ = 7r/2, with the same mass ratio, sin 2 26>^ T = 0.77. 

The value of m V2 /m V3 = 0.05 corresponds to \K\ = 0.63. Since B~ l e l13 = \K\e^ — 2, for 
£ = this gives B = 0.73, or Be — 0.1. In other words, no very great hierarchy is required 
in M R . 

Turning now to case II, we have that 



M," = 



/0 

^ 2 







\ 



e e 
e l + S-V 



i/3 



A 



R 



(29) 



Consequently, for this case 



tan 2 #2 3 



2e/K', 



(30) 



where K' 



1 + B 1 e l/3 . The ratio of eigenvalues of M v gives m U2 /m u . 3 = 
e^y i — \K' — 1\ 2 /\K'\ 2 . If we take, m U2 /m U3 = 0.05, as before, and assume that K' is 
real, we have that K' ^ 0.6. This gives 6% 3 = 12.3°, 6^ ^ 50.5°, and sin 2 26^ ^ 0.96. 
Again, there is good agreement with the SuperKamiokande results. Moreover, since this 
value of K' corresponds to B = —2.5, or Be 2 = —0.05, we see that in this case also no great 
hierarchy is needed in Mr. 

The foregoing discussion is all based on the assumption that the mixing with the first 
family is small, so that one has the small-angle MSW solution to the solar neutrino problem. 
This will certainly be the case if r] = 0. As we will now see, if instead rj = 8 x 10 -6 , as 
needed to have m u = 4.5 MeV, either the small mixing of v e that we have been considering 
or large mixing of v e is possible, depending on the form of Mr. 



B. Mixing of the First Family 

In the previous discussion, we set rj = in which case, no matter what the form of Mr, 
the matrix M v = —^M^N has vanishing first row and column, and the matrix U v that 
diagonalizes M^My has the form of Eq. (26). It is easy to show that the matrix U L which 
diagonalizes LtL has the form 
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cos6»f 9 -sinflfn 



Ur 



sin 6fr 




12 
12 



cos 6h 




/l 

COS $23 
L 

23 



\0 sin 9 



- sin 0f 3 

COS 9^ 3 



(31) 



where sinfl^ = ^t R , t R is defined in Eq. (17), and #23 is given after Eq. (27). Putting these 
together, one has that the total mixing matrix of the neutrinos, U = U L U V , is 



U 



V 



sin 9\ 2 cos O^t — sin 9\ 2 sin 9^ T 
— a C os0f 2 sin0 MT 



I cos 9^ 2 — sin 9^ cos M . 
sin cos 9f 2 cos 



• sin 6^ r 



cos 6> MT 



(32) 



where 9^ T = 9% 3 — #23- This yields the results, already given in Eq. (20), for U e2 and U e3 . 

Now we will consider what happens under what is presumably the more realistic assump- 
tion that r] = 8x 10~ 6 . 

With rj 7^ 0, there are two basic possibilities to consider. One possibility is that Mr has 
a form in which its 12, 21, 13, and 31 elements all vanish or are negligibly small. If such 
is the case, then the previous analysis applies, and the mixing of v e is due entirely to the 
matrix L. The only effect of the parameter 7] in the lepton sector is then to give v\ a mass 
of about 4 x 10~ 7 eV. The second possibility is that Mr does have significant 12, 21 and/or 
13, 31 elements. If this is the case then a strikingly different situation can arise ||, namely 
"bimaximal" mixing |T!|, [ pOj] , 

We will first illustrate what happens with a simple example. Consider the following form 
for M R : 



M 



R 



/ Ae 3 

Ae 3 Be 2 | A R . 
\ 1 



(33) 



We normalize A and B by powers of e simply for later convenience. The mass matrix of 
light neutrinos resulting from this form is 



M v = -N T M^N 



( V 2 B 
1 ^A? 



\ e 2 A 





e 2 
e 



n 1 

e 2 " A 
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1 



Ml 



A 



R 



(34) 



One sees that the 2-3 block has vanishing determinant, so that a rotation in the 2-3 plane 
by an angle #23 — e brings M v to the form 



M' 



( niB_ 
nl 

e A 
_il 

e 2 A 



21 

e A 







Ml 



A 



R 



(35) 



This can be put in a more transparent form by a rotation in the 1-3 plane by an angle 
9" 3 = —r)/(e 2 A). This angle is less than or of order 10~ 4 and thus negligible, practically 
speaking, so 



M" = 




(36) 
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It is clear that the 11 element, being higher order in 77, is likely to be much smaller than the 
12 and 21 elements. The condition for this to be the case is that A/(B — 1) > rj/e 3 = 2 x 1CT 3 . 
If this very weak condition is satisfied, then the form of the matrix manifestly corresponds 
to the situation in which the v e and together form a pseudo-Dirac pair. That in turn 
would mean that the mixing of these two neutrinos is very close to maximal. 

One sees from Eq. (36) that m V3 = M^/Ar, and that the splitting between m Vl and 
m U2 is given by Ara^ = 2(r] 3 /e 5 )((B — l)/A 3 )(Ml I /A R ) 2 . For the vacuum solution to the 
solar neutrino problem, one has Am^ ~ 1CT 10 eV 2 , so that Am^/m^ ~ 3 x 10~ 8 = 
2{rf/e b ){B — 1) /A 3 . This gives A(B — I) 1 / 3 0.06. Thus no great hierarchy is required in 
M R to get the vacuum oscillation solution. The reason for this is that in this scheme the 
smallness of Am^ is due to the extreme smallness of the parameter 77, which is equal to the 
ratio m u /m t . 

It is easy to see from what has already been said that the matrix U v needed to diagonalize 
MlM v is of the form 

/ 1 \ / l/y/2 1/V2 0\ 

U v 2* cos 0^ sin 0^ -l/y/2 l/y/2 . (37) 
\0 -sin ^ 3 cos0£ 3 / V 1/ 

where we have neglected the tiny rotation 9^. The matrix Ul is already given in Eq. (31), 
so that the full neutrino mixing matrix can be written 

/ l/y/2 l/y/2 0\ 
U = U S ma- -l/y/2 y/2 , (38) 

V 01/ 

where Usma is given in Eq. (32), and is just the form that results in the small mixing angle 
(SMA) MSW case of this model. In other words, the net result of the large mxing of the 
first family produced by the A entries in Eq. (34) is simply to multiply the (SMA) MSW 
form of U on the right by a rotation of 7r/4 in the 1-2 plane. Consequently, the predictions 
for C/ M 3 and U eS are essentially unaffected. However, U e2 becomes 1 / y/2 instead of the value 
given in Eq. (20). 

The interesting lesson is that "bimaximal mixing" is easy to achieve if the large mixing 
of and v T comes from the charged lepton sector, i.e., from diagonalizing L, while the large 
mixing of v e comes from diagonalizing M v . 

The simple form given in Eq. (33) gives 0£ 3 ^ e = 8°, and thus 9^ T = 9^ - 6^ = 55°, 
corresponding to sin 2 29^ T = 0.88. Somewhat larger values of sin 2 29^ T can arise from a more 
general form Mr. Consider, for example, 



M R = Ae 3 Be 2 A R . (39) 




Then 

/ vIb V-BC v(BC-A) 



"BU *{BU-A)\ M2 
M v = - z BC e 2 A 2 eA(A + C) \{A 2 + BC 2 )- 1 -^ (40) 



\Z(BC-A) eA(A + C) (A + C) 



2 
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Here, as in Eq. (34), the 2-3 block has vanishing determinant. The crucial difference is that 
the diagonalization of this matrix involves a rotation in the 2-3 plane by an angle 9% 3 = t-Ajj- 
With C = —\A, for instance, 6^ T comes out very close to 45°. Otherwise, this case is quite 
similar to that of Eq. (33). 



VI. DETAILS OF A SPECIFIC MODEL 

In the previous Sections we have presented the construction of our SO (10) minimal Higgs 
model in the framework of effective 5*0(10) and 577(5) operators. We now show that one 
can construct a technically-natural realization of this scheme by introducing sets of Higgs 
and matter superfields with a well-defined family symmetry. We first address the Higgs 
sector. 



A. Higgs Sector with U{1) x Zi x Z2 Family Symmetry 

The doublet-triplet splitting problem in SU(5), and therefore 50(10), arises because 
the colored Higgs in the 5 — 5 pairs of each 10# must be made superheavy at the GUT 
scale, while just one pair of Higgs doublets should remain massless there and be free to 
develop VEV's at the electroweak scale. This problem has been addressed and solved in || 
by the introduction of just one 45 adjoint Higgs field with its VEV pointing in the B — L 
direction, together with two pairs of 16 + 16 spinor Higgs fields, two 10 Higgs in the vector 
representation plus several Higgs singlets. We shall briefly summarize the solution, but first 
we note that it is necessary to introduce several more Higgs fields in the vector and singlet 
representations in order to generate the Yukawa structure for the fermion mass matrices 
presented earlier. 

The authors of || found that the Higgs superpotential required to solve the doublet- 
triplet splitting problem could be neatly obtained from their list of Higgs fields by introducing 
a global family symmetry group of the type U(l) x Z2 x Z2, which can arise in a natural 
fashion from string theory. With this in mind, we now list in Table I all the Higgs fields to 
be considered together with their family charge assignments. 
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Higgs Fields Needed to Solve the 2-3 Problem: 



45 B 
16 
16 
10 
1: 



A(oy 

cay 



x(oy 



pY 



++ 



, c>(-l-p)-+ 

T 2 (-l)+- 

P(p) + -, Z 2 (p) 
Additional Higgs Fields for the Mass Matrices: 
10: 



1: 



T (l+p) + -, 72(l + 2p) + - 
T (-3+p)-+,^(-l-3p)- 
F(2)-+, Y'(2)++, 5(2 - 2p)" 
Vm(4 + 2 P )++ 



S"(2 - 3p)" 



Table I. Higgs superfields in the proposed model. 

As noted in the table, in order to complete the construction of the Dirac mass matrices, 
four more vector Higgs fields and four additional Higgs singlets are needed, while one Higgs 
singlet is introduced to generate the right-handed Majorana neutrino mass matrix. 

It is then possible to write down explicitly the full Higgs superpotential from the Higgs 
50(10) and family assignments, where we have written it as the sum of five terms: 



Wffiggs 

W A 
W CA 

W 2 / 3 

Wr 



W A + W CA + W 2/3 + W Hd + Wr 

trA A /M + M A trA 2 

X(CC) 2 /M 2 C + F{X) 

+C'(PA/M 1 + Z X )C + C(PA/M 2 + Z 2 )C 

T 1 AT 2 + Y'T 2 



(41) 



T^CCY'jM + T CC + T (T S + T^S') 
TqTqVm 



The Higgs singlets are all assumed to develop VEV's at the GUT scale. We can then 
determine the fate of the other Higgs fields from the F-flat and D-flat conditions. In par- 
ticular, Wa fixes (A) through the Fa = condition where one solution is (A) oc B — L, the 
Dimopoulos-Wilczek solution . Wqa gives a GUT-scale VEV to C and C by the Fx = 
condition and also couples the adjoint A to the spinors C, C, C and C without destabiliz- 
ing the Dimopoulos-Wilczek solution or giving Goldstone modes, as shown in ||. W 2 /% gives 
the doublet-triplet splitting by the Dimopoulos-Wilczek mechanism Jl2|, [|]. Wh d mixes 



the (1, 2, —1/2) doublet in T\ with those in C (by Fq = 0), and in T and Tq (by Ff o = 0). 



B. Yukawa Sector 

We now turn to the Yukawa sector and specify the matter fields and their U(l) x Z 2 x Z 2 
charge assignments which will complete the realization of the specific model in question. 
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For this purpose, we require three spinor fields 16«, one for each light family, two vector-like 
pairs of 16 — 16 spinors which can get superheavy, a pair of superheavy 10 fields in the 
vector representation, and three pairs of superheavy 1 — l c fermion singlets. The complete 
listing is given in Table II. 

l 6l (-|-2p)+- 16 2 (-i + p)++ 16 3 (-±) ++ 

16H i6'Hr + 

16(f )+- T6'(-| + 2p) + ~ 

10i(-l-p)- + 10 2 (-1+J9)++ 

l x (2 + 2p)+- l a (2-p)++ 1 3 (2) ++ 

i;(-2-2p)+- l c 2 (-2)+- m- 2 -p)++ 

Table II. Matter superfields in the proposed model. 

In terms of these fermion fields and the Higgs fields previously introduced, one can then 
spell out all the terms in the Yukawa superpotential which follow from their SO (10) and 
U(l) x Z 2 x Z 2 assignments: 

W Yukawa = 16 3 • 16 3 • Tx + 16 2 • 16 • Ti + 16' • T x 
+16 3 • 16i • Tq + 16 2 • 16i • T + 16 3 • 16 ■ A 
+16i • Ttf ■ Y' + 16 • 16 • P + 16' • IB' ■ S 

+16 3 ■ 10 2 _C" + 16 2 ■ 10j_ C + 10i • 10 2 ■ Y (42) 
+16 3 ■ 1 3 ■ C + 16 2 ■ 1 2 • C + 16i • li • C 
+1 3 -1 C 3 -Z + 1 2 -1 C 2 -P + 1 1 -1 C 1 -X 
+11 ■ 11 ■ V M + 11-11- V M 

where the coupling parameters have been suppressed. To obtain the GUT scale structure 
for the fermion mass matrix elements, all but the three chiral spinor fields in the first line of 



Table II. will be integrated out to yield Froggatt-Nielsen diagrams [^] of the type pictured 
earlier. Note that the right-handed Majorana matrix elements will all be generated through 
the Majorana couplings of the Vm Higgs field with conjugate singlet fermions given in the 
last two terms of Eq. (42). 

In order to present a clearer description of how the GUT scale mass matrices are deter- 
mined from the Yukawa and Higgs superpotentials, we shall illustrate the procedure for the 
up quark mass matrix U. The three massless color-triplet quark states each with charge 2/3 
are obtained as linear combinations of all such color and charge states within the fermion 
supermultiplets given in (41). In particular, the basis for the left-handed states ul and 
left-handed conjugate states u c L can be ordered as follows: 

&u L = ||3, 2, | >io(i6i)j |3> 2, \ >io(i6 2 )> |3, 2, | >i (i6 3 ), |3, 2, | >io(ie) ■> ^ 
|3, 2, | >io(i6')> |3, 1, | >io(l6)> |3, 1, | >To(T6')} 

Bu c L = {|3, 1, — | >10(16i)j |3, 1, — | >10(16 2 )> |3, 1, — | >10(16 3 )) |3, 1, — | > 10(16) ) 



|3, 1, — | >io(i6')> |3, 2, — | >to(T6)' |3, 2, — | > io(i6')} 



(44) 
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where the states are labeled by their representations and hypercharge according to 
\SU(3) C , SU(2) L , Y) 

Sf/(5)(5O(10))- 

We then form the Yukawa contribution u c ^C~ x D u u L by using the above bases and the 
superpotentials to obtain the matrix 



D„ 




















y'\ 











t 2 

























a 








t 2 











P 

















t' 





s" 








a 


p 











\y' 











s" 





o) 



where we have introduced the following shorthand notation: 



a 



A163I63T1 (Ti) 

^16 3 T6a(^)? 
16'T6'5 X 



h — Ai6 2 i6Ti(Ti), t' — XlQ'WT! (Ti) , 



p 
y' 



1616P 



I61I6 Y 



(P), 



(45) 



We can then determine from this matrix the three pairs of zero-mass eigenstates at the GUT 
scale where the electroweak VEV of Ti vanishes: 



Uil) = 


110(160) 


u 2 l) = 


10(16 2 )> 


U3L) = 


|10(16 3 )> 


«?L> = 


110(160) 


u °2l) = 


10(16 2 )) 


U h) = 


"|10(16 S )> 



(46) 



and where the states are now simply labeled by their SU(5) and 5*0(10) representations. 

Finally, the Dirac matrix U for the three light quark states u, c, t is obtained by 
bracketing the electroweak contributions by the appropriate u\ L state on the left and the 
Uj L state on the right. The result obtained for U has exactly the form found earlier in Eq. 
(15) from the previous effective operator approach, with the identifications: 



Mu= (* 3 ) 5 (io) 
eM v = |3(a 9 /p)(i 2 ) 5 (io)| 
V M V = (y'/s"f(t% m 



(47) 



Here the subscript on a q signifies a factor of 1/3 arising from the B — L VEV of the A in the 
adjoint representation, while the subscripts on the t terms specify the appropriate doublet 
VEV in the 10 for Ti. We have neglected the state normalization factors in (47) but will 
later argue that they can all be taken to be approximately unity. 

The Dirac matrices, D, N, L are constructed in a similar fashion. In the case of D 
and L, the bases corresponding to Eqs. (43) and (44) are enlarged by two states lying in 
the 10i and 10 2 representations of 5*0(10). For N, on the other hand, in addition to the 
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two above states, one must add the singlet fermion contributions from the representations 
1^, k — 1, 2, 3 for B UL and 1%, /c = 1, 2, 3 for B v ^. We list the zero-mass state vectors for D, 
N and L in analogy to Eqs. (46): 



\d lL ) = 
\d 2L ) = 
\d 3L ) = 
\d\ L ) = 
\d c 2L ) = 
\dh) = 



|10(16i))-^|10(16'))J /^l + y' 2 /s" 2 
10(16 2 )> 

|10(16 3 ))-f|10(16))] /y/l + aVp 2 



15(160) - JL| 5 (16')> /Jl 



yl2/ s >,2 



|5(16 2 )) - ||5(10 2 ))] / 

|5(16 3 )> - § |5(16)>] //T + ^y 



(48) 



n 2 L) 

) 
) 



n 



1L 



n 



2L 



n 



3L/ 



^|5(i6'))J/Vi + y'7s" 2 

5(102)>1 /v/l + c 2 /?/ 2 



15(160) 
15(162)) - J 

|5(163)) -5|5(16))l/yi + aVp 



1(160) " £|1(16')) " till) /V 1 + ^/s" 2 + c 2 /x 2 

i(i60)-^|i^)lA 



(49) 
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l(16 3 ))-f|l(16))-f|l 



1 + a 2 /p 2 + cl/z 2 



hi) 

hi) 



<§L> 



|5(160) - JL|5(16')) /\/l+y' 2 /s" 2 
|5(16 2 )> 

|5(16s)) 



f 15(102)) /yi + cVy 2 



5(16)) /Jl + a*/p 



|10(160) - £|10(16'))] /yjl + y^/s 
10(16 2 )) 



(50) 



|10(16 3 ))-S|10(16)) /Jl + ayp 



In the above we have introduced, in analogy with Eqs. 
notation: 



(45), the additional shorthand 



c — Al6 2 10iC\C); 

x = A lll c X (X), 



Xi 3 ic z(Z), 



Ci 

y 

P22 



A 



(C),i = 1,2,3, 



A iOiio 2 y 

'Ws-p(-P) 



(51) 



The Dirac matrices D, AT and L are found by forming matrix elements of the electroweak 
symmetry breaking VEV's with the appropriate basis vectors. Again these matrices have 
exactly the structures given in Eqs. (15), provided the state normalization factors are 
approximated by unity, i.e., we assume that the zero-mass states have their large components 
in the chiral representations 16i, 16 2 and 16 3 and that all the other components are small. 
We shall return to this point in the next Section. In the meantime we make the identifications 
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M D 
eM D 
a Mr, 
AM D 
6'e^M D 



(£3)5(10) 

|3(a g /p)(t 2 )s(io)| 
-(c/y)(c')- 5(w) 

toh/s 



(52) 



in terms of the notation given in Eqs. (45) and (51) and the following: 



to 
to 
s 



^16il6 2 T ? 

A CC'T (C)(C>), 



t'o= A 
c'= X 



16il6 3 T^ 
I63IO2C 



(C), 



(53) 



A 



(S>) 



The phase appearing in the 5' term can be understood to arise from a phase in the VEV 
for S'. The structures of the Dirac matrix elements given in Eqs. (15), (47) and (52) can 
be understood in terms of the simple Froggatt-Nielsen diagrams of Fig. 1 and 2, with the 
Higgs fields labeled as in Table I. 

Turning to the right-handed Majorana mass matrix, we use the zero mass left-handed 
conjugate states that were obtained implicitly above for the Dirac matrix N to form the basis 
for Mr. The right-handed Majorana matrix is then obtained by bracketing the Majorana 
Higgs Vm with the appropriate zero mass conjugate neutrino states in (49). We obtain 



M, 



R 



I 

V 



Ae 3 





0\ 


1/ 



A 



R 



(54) 



where 



Mo 



M 3 = A R 
-M l = Ae 3 A R 



Xic 1 c Vm (V m )(c 3 /z) 2 , 

^imV M ( V M)(c 1 /x)(c 2 /p22) 



(55) 



The lighter two right-handed Majorana masses are degenerate and have opposite CP-parity. 
Note that the whole right-handed Majorana matrix has been generated in this simple model 
by one Majorana VEV coupling superheavy conjugate fermion singlets as shown in Fig. 4. 

We conclude this Section with a summary of the GUT scale predictions derived from the 
Dirac and Majorana mass matrices with the particular parameters appropriate for the model 
in question. For convenience we give the whole set equations which are the counterpart of 
Eqs. (19). 
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(a 2 + 1; 

¥ 2 ■ [i - 

If _2_ . 



m°/m° f L f fl • [1 

'A 

ml/ml ^ 



-y 2 Mu/M D , m°Jm° t rj, 
§e 2 ], ml/m^l-l^ 



\1 



1 l-CT 2 -g-£/3 

3 6 ct(<t 2 +1) 



i e (T (^2 4 f 2 i) (<L + t 2 R ) 



3 o-2+l 6 ' 
5 (*£ + *£)]> 



u° ^ 

y us 



u° e2 = 



ttO ~ 

^e3 — 



e ^ 2 TI ' t 1 + e CT (a2+i) £ + 18 + *Jl)]) 

3 e ?2+l ' t 1 + 3 e a(<r 2 +l)]' 



i j 7 



hi 



1 



1 - 



vV 2 


+ 1 


1 


1 






+ 3Vc 


1 


( 


3vV 2 +l ^ 



a — e 



_2 



a 2 + l 



3a " L 



o" - e)t R - 



At 2 



(56) 



To the quark equations we have added the ratio m°/m° which involves the ratio of (5(Ti)) 
to (5(Ti)), i.e., Mjj/M d , as well as giving the leptonic mass ratios and mixings specific to 
the model in question. 



C. Numerical Evaluation of Matrix Parameters 

We have elaborated above how the simple explicit model proposed gives precisely the 
structure for the Dirac mass matrices that was obtained from the effective operator approach. 
We now show that the entries are also numerically in the range to fit the quark and lepton 
mass and mixing data. 

In order to compare the GUT scale predictions in Eq. (56) with the low scale data, the 
GUT scale values are first evolved from A G = 2 x 10 16 GeV down to the SUSY scale which 
is taken to be A SUSY = m t (m t ) = 165 GeV by use of the 2-loop MSSM beta functions. For 
this purpose, the mass ratios at the two scales are related by the rji/j running factors, while 
the quark mixing elements are scaled by the rjij factor according to 

( m i \ 

V m j ) SUSY 
(^j)sUSY 



In. 



'i/ji 



VfJVij, (u) = (ub), (c6), (td), (ts) 



(57) 
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The remaining evolutions to the bottom and charm quark or tau lepton running mass scales, 
or to the 1 GeV scale in the case of the light quarks and leptons, is carried out with the 
3-loop QCD and 1-loop QED renormalization group equations. Here the running factors are 
r/i with the mass ratios scaled according to 

mi(mi) = (mi) SVS Y /Vi( m t) (58) 

or similarly, with the running mass scale TOj replaced by 1 GeV. With tan ft = 5 used 
for the numerical evaluations for reasons that will become apparently shortly, a s (Mz) — 
0.118, a(Mz) = 1/127.9 and sin 2 9\y = 0.2315, the running factors are given by 

Vu/t = Vc/t = 0.6927, i] d/b = 7] s/b = 0.8844 

■np/r = 0.9988, Vb/t = 0.5094 

Vub = Vcb = Vtd = rjts = 0.8853 



r) u {mt) = 0.4235, r] c (m t ) = 0.4733, r] t (m t ) = 1.0000 
Vd (m t ) = 0.4262, r] s (m t ) = 0.4262, r] t (m t ) = 0.6540 
rj e (mt) = 0.9816, r^(m t ) = 0.4816, r) r {m t ) = 0.9836 



(59) 



Finally, finite corrections must be applied to m s , m b and the evolved quark mixing matrix el- 
ements which arise from gluino and chargino loops. The correction factors are conventionally 
written as (1 + A s ), (1 + A b ) and (1 + A cb ) where we have used 

A s = -0.20, A b = -0.15, A cb = -0.05 (60) 

as explained below. 

Using the quantities m t {m t ) = 165 GeV, m T = 1.777 GeV, = 105.7 MeV, m e = 



0.511 MeV, m u = 4.5 MeV, V us = 0.220, V cb = 0.0395, and 5 CP = 64° to determine the 
input parameters, one obtains for them Mu ~ 113 GeV, Md — 1 GeV, a = 1.780, e = 
0.145, t L = 0.236, t R = 0.205, 9 = 34° (corresponding to 5 = 0.0086, 5' = 0.0079, = 54°), 
and rj = 8 x 10~ 6 . With these inputs the remaining quark masses and mixings are obtained, 



to be compared with the experimental values [22] in parentheses: 



m c (m c ) = 1.23 GeV (1.27 ±0.1 GeV) 

m b {m b ) = 4.25 GeV (4.26 ± 0.11 GeV) 

m s {l GeV) = 148 MeV (175 ± 50 MeV) (61) 

m d (l GeV) = 7.9 MeV (8.9 ± 2.6 MeV) 

\V ub /V cb \ = 0.080 (0.090 ± 0.008) 

where the finite SUSY loop corrections for m b , m s and V cb have been rescaled to give 
m b (m b ) ~ 4.25 GeV for tan/3 = 5. Had we chosen Sep = 70° as input, on the other hand, 
we would find instead (V^/V^l = 0.085. With the numerical values in (61) we find for p, fj 
and the a, [3 and 7 angles of the unitarity triangle pictured in Fig. 3 

p = 0.143, 77 = 0.305, a = 96°, (3 = 20°, 7 = 64° (62) 

The upper vertex of the triangle appears to be circled precisely in the allowed experimental 
region. 



28 



Additional predictions follow for the neutrino sector. The effective light neutrino mass 
matrix of Eq. (34) or (36) with B = leads to bimaximal mixing with a large angle solution 



for atmospheric neutrino oscillations || and the "just-so" vacuum solution [B5| involving 
two pseudo-Dirac neutrinos, if we set A R = 2.4 x 10 14 GeV and A = 0.05. We then find 

m 3 = 54.3 meV, m 2 = 59.6 /xeV, mi = 56.5 fieV 
M 3 = 2.4 x 10 14 GeV, M 2 = M 1 = 3.66 x 10 10 GeV 

> 

(63) 





U e2 = 0.733, U e3 = 


0.047, C/ M3 = 


-0.818 


5 S'CP 


= -0.2° 


Am 2 3 


= 3.0 x lO" 3 eV 2 , 


sin 2 20 atm 






2 = 0.89 


Am 2 12 


= 3.6 x 10- 10 eV 2 , 


sin 2 2# solar 


= Mu el 


\ 2 We2 


2 = 0.99 


Am 2 3 


= 3.0 x lO- 3 eV 2 , 


sin 2 26> rea c 


= 4|f/ e3 


2 (1- 


\U e3 \) 2 = 0.009 



The effective scale of the right-handed Majorana mass contribution occurs two orders of 
magnitude lower than the SUSY GUT scale of Ac = 1.2 x 10 16 GeV. The effective two- 
component reactor mixing angle given above should be observable at a future neutrino 



factory, whereas the present limit from the CHOOZ experiment jzo] is approximately 0.1 for 
the above Am 2 3 . In principle, the parameter A appearing in Mr can also be complex, but 
we find that in no case does the leptonic CP-violating phase, 5' CP exceed 10° in magnitude. 
Hence the model predicts leptonic CP-violation will be unobservable. 

The vacuum solar solution depends critically on the appearance of the parameter 77 in 
the matrix N, corresponding to the non-zero r\ entry in U which gives the up quark a mass 
at the GUT scale. Should we set T] = 0, only the small-angle MSW solution [|7| would be 
obtained for the solar neutrino oscillations. The large angle MSW solution is disfavored by 
the larger hierarchy, i.e., very small A value, required in Mr. 

Finally we address the issue that the state normalization factors were all replaced by 
unity in Eqs. (47) and (52) for the various matrix parameters. This is a good approxi- 
mation provided the three fermion spinor states |16i), 1 162) , 1 163) provide the dominant 
contributions to the zero-mass quark and lepton states at the GUT scale. In particular, the 
following ratios must be much less than unity: 

(a/p) 2 , {y'/s")\ {c/y)\ (cjxf, (c 2 /p 22 ) 2 , {c 3 /zf « 1 (64) 

Let us assume for simplicity that the electroweak couplings of (Ti) in t 3 , t 2 and t' in Eq. 
(45) and of (C) in d of (53) are identical. Then with e = \3(a q /p) \ = \a e /p\ = 0.14, we find 
the condition (a/p) 2 ~ 0.02 1 holds. To obtain an up quark mass m u (lGeV) ~ 4.5 MeV, 
we need 77 = (y'/s") 2 ~ 8 x 10~ 6 at the GUT scale, which easily satisfies (64). 

Requiring that (c/y) 2 1 and with the result from Eqs. (52) that 



a 



c (5(C)) 

y (5(2i)) 



1.8 (65) 



leads us to the results that 

(66) 



tan/3 ~ \J o 2 + lm^(cos'y) /m° mf/mjj 
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in terms of the T\—C mixing angle, 7, in Eq. (11). With c/y = 0.1, for example, tan7 ~ 18 
which implies tan/3 ~ 6, a very reasonable mid-range value allowed by experiment. For this 
reason, we have chosen to illustrate the numerical results above with tan /3 — 5. 

The remaining ratios in Eq. (64) can also be satisfied. For comparable q's, Ae 3 ~ 

1.4x 10 -4 obtained from Eq. (55) requires that (Z) / (X)(P) ~ 0.01. This ratio is consistent 
with the VEV's needed in the Higgs superpotential of Eq. (41) in order to solve the doublet- 
triplet splitting problem. 

Turning now to the parameters 5 and 5', we note that the near equality of their mag- 
nitudes leads to the ratio S/\S'\ = s/\s'\ — 1. Moreover, if we assume y ~ y', we obtain 
the estimate 5 ~ cy'/(ys) tan 7 ~5x 10~ 3 with the numbers obtained earlier, whereas the 
actual value needed is 5 ~ 0.008. 

Thus we have found that not only are the desired forms of the Dirac (and Majorana) 
matrices generated by the model of this Section, but that the numerical values required for 
the matrix parameters are also quite reasonable. 

VII. SUMMARY 

Both the largeness of the atmospheric neutrino mixing and the smallness of the quark 
mixing V c b can be elegantly accounted for by the idea that the charged lepton mass matrix L 
is highly asymmetric or "lopsided" and that the down-quark mass matrix D is related to the 
transpose of L by an 577(5) symmetry. This idea was discovered independently by several 
groups and has since been used in numerous models of fermion masses. Remarkably, exactly 
such mass matrices emerged in our work from quite other considerations than neutrino 
masses and mixings, specifically from an attempt to construct the simplest possible realistic 
5O(10) model. 

Advances have been made in recent years in simplifying the Higgs structure of SUSY 
5O(10) models. If one assumes the minimal set of Higgs fields that can break 5O(10) down 
to the standard model group, the possibilities for Yukawa terms for the quarks and leptons 
become significantly restricted. It turns out that there is what seems to be a uniquely simple 
set of 5O(10) Yukawa terms that gives realistic masses and mixings. This set consists of 
only six effective Yukawa terms (five if m u = 0) which satisfactorily fits all nine masses 
of the quarks and charged leptons as well as the four CKM parameters. In addition, large 
— v T mixing emerges automatically. Moreover, in this uniquely simple model, the simplest 
possibilities for the Majorana mass matrix M R of the right-handed neutrinos lead either to 
small angle MSW values for the solar neutrino mixing or to vacuum oscillation values. In this 
paper we have studied in detail the consequences of different forms of Mr for the neutrino 
mixing angles and mass ratios. 

In the published literature no more predictive and economical a model for quark and 
lepton masses than the one discussed here exists that is also consistent with present knowl- 
edge. It is striking that in this model a single term and a single parameter (which we call a) 
accounts for no less than four puzzling aspects of the light fermion spectrum: the largeness 
of U^, the smallness of V c b, the smallness of m c /m t compared to m s /mb, and the Georgi- 
Jarlskog factor of three between and m s at the GUT scale. It should be emphasized that, 
while many satisfactory neutrino mixing ideas and also many interesting ideas for explaining 
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the pattern of quark and charged lepton masses have been proposed, very few models exist 
which not only give a satisfactory account of neutrino phenomenology but are at the same 
time highly predictive. 

We have shown that the model defined by the existence of these five (or six) effective 
Yukawa terms can be realized in a complete and specific renormalizable SUSY SO (10) 
model that is technically natural. We have presented the details of such a model, including 
all the Higgs and quark and lepton superfields, the abelian flavor symmetries, and the 
transformation properties of the fields under these symmetries. Finally, we have done a 
quantitative comparison of the predictions of the model to experiment. 

In the future this model will be rigorously testable in several ways. The most important 
are (1) a relation between the real and imaginary parts of V u b including a precise test of the 
angles of the unitarity triangle; (2) a prediction for U e 2, which in the small angle MSW case 
gives a sharp relation between the solar and atmospheric angles; and (3) a definite prediction 
for U e3 . 
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FIG. 1. Diagrams that generate the elements in the quark and lepton mass matrices shown 
in Eqs. (10) with the Higgs labeling corresponding to that appearing in Table I of Section VI. 
(a) The 33 elements denoted "1". (b) The 23 and 32 elements denoted "e". Note that because 
of the appearance of the VEV of the adjoint Higgs field 45h = A, they are proportional to the 
50(10) generator B — L. (c) The asymmetric elements denoted "cr" arise from this diagram. That 
they do not contribute to the up quark masses, and contribute asymmetrically to the down quark 
and charged lepton mass matrices are consequences of the fact that the 5O(10) 10's, i.e., 10i and 
10 2 , contain 5 but not 10 of 5*7(5). 
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FIG. 2. Diagrams that generate the masses of the first family of quarks and leptons. See Eq. 
(14). (a) The 11 element called (b) The 12 and 21 elements called " <5" . Because the 10 H 

couples to the symmetric product of I61I62, S appears symmetrically in the mass matrices, (c) 
The 13 and 31 elements called "(5"', which also appear symmetrically. 
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FIG. 3. The unitarity triangle for V u dV* b + V c dV* b + VtdVtb = is displayed along with the 
experimental constraints on V U dV* b , which is the upper vertex in the triangle. The constraints 
following from \V ub \, B-mixing and e extractions from experimental data are shown in the lightly 
shaded regions. The experimentally allowed region is indicated by the darkly shaded overlap. The 
model predicts that V u dV* b will lie on the dashed circle; cf. Eq. (21). The particular point on this 
circle used to draw the triangle shown is obtained from the numbers given in Section VI; cf. Eq. 
(62). 
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FIG. 4. Diagrams that generate the 33, 12, and 21 elements of the Majorana mass matrix Mr 
of the superheavy right-handed neutrinos. 
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